Enhanced laboratory sensitivity to variation of the fine-structure constant using 

highly-charged ions 
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We study atomic systems that are in the frequency range of optical atomic clocks and have 
enhanced sensitivity to potential time- variation of the fine structure constant ct. The high sensitivity 
is due to coherent contributions from three factors: high nuclear charge Z, high ionization degree, 
and significant differences in the configuration composition of the states involved. Configuration 
crossing keeps the frequencies in the optical range despite the large ionization energies. We discuss 
a few promising examples that have the largest a-sensitivities seen in atomic systems. 

PACS numbers: 06.30.Ft,37.10.Ty,31.15.am,32.30.Jc 



I. INTRODUCTION 

Theories that seek to unify gravity with the other fun- 
damental interactions suggest that temporal variation of 
fundamental constants is a possibility, or even a neces- 
sity, in an expanding Universe (see, e.g. [l|). Hints from 
quasar absorption spectra that the fine-structure con- 
stant, a — e^ /fic, may have been different in the dis- 
tant past have been reported [2, Q but have not been 
confirmed by other groups working on a different tele- 
scope Q. More recently, the methodology of [4| was 
questioned Q and a reanalysis of the same data increased 
the reported error bars by a factor of six [6| . The obser- 
vational status is therefore still unclear. 

Atomic clocks provide a complementary method to 
search for temporal variation of fundamental constants 
in terrestrial laboratories [3, @] . Two different clocks are 
compared over the course of several years and any dis- 
crepancies are interpreted in terms of variation of a. For 
such a scheme it is critical that the two clocks have dif- 
ferent sensitivities to any possible variation of a. The 
best current laboratory limit, a/a = (—1.6 ± 2.3) x 
10"^'' year~^, comes from comparison of Hg+ and A1+ 
optical clocks over the course of a year [9| . Here the Hg+ 
transition is strongly dependent on a, while the A1+ clock 
is practically insensitive 0, [l3 ■ 

The laboratory clock limits on a may be improved by 
finding systems with enhanced sensitivity to a- variation, 
usually denoted q (see Section |ll|. Candidates under 
consideration include the thorium "nuclear clock" [ll|, 
which would utilise the ~ 7.5 eV nuclear transition in 
229rp-|^ 12]. Among other benefits, such a clock could have 
many orders-of-magnitude larger g- values than the opti- 
cal Hg+ clock transition [13, 14]. Optical atomic transi- 
tions with high sensitivity can be found in Yb+ [15| and 
Th3+ 16]. 

A different kind of sensitivity comes from transitions 
that have high relative sensitivity to a. For example, 
in the dysprosium atom there are two different transi- 
tions that are "accidentally" ne arly degenerate and have 
g- values with different signs [3, [13, Ea] • Because the lev- 
els are so close together the relative sensitivity, defined 



hy K = 2Ag'/w where uj is the frequency of transitions 
between the two levels, is extremely high (~ 10^). The 
first experiment to utilise this transition gave fairly tight 
limits on a- variation [19|, but the full enhancement was 
not realised because one of the levels is very broad. 

In this paper we show that both kinds of sensitivity 
(large q and large K) can be realized in highly-ionised 
atomic systems. While atomic spectroscopy in electron 
beam ion traps is currently not competitive with optical 
frequency standards (see, e.g., [20, [2l[ and review [22]) 
the technology continues to improve, and with the en- 
hancements in sensitivity reported here, highly-charged 
ions may prove to be a good system for detecting vari- 
ation of a. In Section HI] we show, using the Ag isoelec- 
tronic sequence as an example, why high q- values can oc- 
cur in highly-charged ions, and how the tendency of such 
systems towards large transition frequencies can be over- 
come. In Section IIVI we show the results of our atomic 
calculations (Section IIII[) applied to some of the most 
promising ions in the sequence. In addition we identify a 
two- valence-electron ion, Sm^^+ , which has optical tran- 
sitions that are the most sensitive to potential variation 
of a ever found. 



II. THEORY 

Using a simple analytical estimate of the relativistic 
effects in transition frequencies, we can see that more 
highly charged ions have higher sensitivity to a varia- 
tion. Note that a ratio of frequencies (the quantity that 
is actually measured) does not depend on the units one 
uses. In this paper we use atomic units e = rrie = Ti = 1 
unless otherwise stated; in these units the atomic unit of 
energy is constant. Consider the relativistic corrections 
to the central-field Schrodinger equation for a valence 
electron, derived in the Pauli theory (see, e.g. j23| ) 
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which is accurate to order {v/c)"^. Here R{r) and Eq 
are the non-relativistic radial wavefunction and energy, 
respectively, V{r) is the potential, and X = jij + l)—l{l + 
1) — s(s + 1). Near the origin V — Z/r (it is unscreened 
by core electrons) and the integrals ([T]) converge as 1/r^. 
The radial s-wave wavefunction near the origin is given 
in the semiclassical approximation by (see, e.g. [2J|) 
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where Za is the effective charge that an external electron 
"sees" and ly is the effective principal quantum number, 

defined by Eg — —jfi- For a single valence electron 
above closed shells, Za = Zi + 1 where Zi is the ion 
charge, while for hydrogen-like ions Za = Z and v — n, 
the principal quantum number. For higher waves Rf{r) 
is proportional to the same parameter Z^Zjv^ . 

Noting that the integrals converge over a distance 



< 



3/(2Z) and neglecting the small contribution of the 



£^0 term in (P (for v > Zl/Z'^) we obtain 
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where /„ is the ionization energy of the orbital. Actu- 
ally, if we neglect the Eq term, we can immediately see 
from ([T]) that the ratio of the relativistic corrections in 
our "distorted wave" case to the relativistic corrections 
in the Coulomb case is proportional to the ratio of the 
squared wavefunctions near the origin: 
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This again leads to Eqs. (|31) and (g]). 

In principle we could restore the Eq term in ([T]). How- 
ever many-body corrections are much more important 
than small improvements to Eq. ([3]). As discussed in 
[lO| , for single- valence electrons the relativistic energy 
shift can be approximately described by the equation 
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where C{Z,jJ) depends on the atom and partial wave, 
but does not depend strongly on the principal quantum 
number. 

In practice a more sophisticated numerical treatment 
is needed. We characterise the dependence of transition 
frequencies on small changes in a by the parameter q, 
defined by the formula 



U}Q + qx 
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where ujq is the energy at the present day value of the 
fine-structure constant ao and 
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FIG. 1. Ratio —q/{Z^a^In) for the 5s level along the Ag iso- 
electronic sequence, calculated using the Dirac-Fock theory. 
It is seen that as Z increases the ratio tends to a constant 
(which is close to 1/5) in accordance with ([3]). 



Our definition of q refiects the fact that the relativistic 
shift scales as a^ . Calculations of cu and q are discussed 
in Section IhII 

When our formula (JH) is correct (when many-body ef- 
fects are not too large), q w A„/a^. In this case, the 
ratio —q/{Z'^a'^In) is a constant. In Fig. [T] we present 
a Dirac-Fock calculation of this ratio for different ions 
along the Ag isoelectronic sequence: Ag, Cd^, In^+, etc. 
It is seen that as Z increases the ratio indeed tends to a 
constant (which is close to 1/5). 

We have shown that sensitivity to a- variation increases 
with ion charge as Z^ = {Zi + 1)^. Unfortunately, the 
interval between different energy levels in an ion also in- 
creases as ~ Za, which can quickly take the transition 
frequency out of the range of lasers as Za increases. How- 
ever, the phenomena of Coulomb degeneracy and config- 
uration crossing can be used to combat this tendency. In 
a neutral atom, an electron orbital with a larger angular 
momentum is significantly higher than one with smaller 
angular momentum but with the same principal quan- 
tum number n. On the other hand, in the hydrogen-like 
limit orbitals with different angular momentum but the 
same principal quantum number are nearly degenerate. 
Therefore somewhere in between there can be a crossing 
point where two levels with different angular momentum 
and principal quantum number can come close together: 
in such cases the excitation energy may be within laser 
range. 

Consider again our example of the Ag isoelectronic se- 
quence. Neutral Ag {Z — 47) has a single valence elec- 
tron above closed shells. The ground state has the va- 
lence electron in the 5s orbital, while the 4/ orbital forms 
an excited level. In Fig. [3] we present calculated Dirac- 
Fock ionisation energies of the Ag isoelectronic sequence. 
One can see that as Z is increased, there is a crossing 
point where the 4/ level becomes the ground state. At 
this point, around Z = 61, even though both levels have 
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TABLE I. Energies and sensitivity coefficients (g) for isoelec- 
tronic sequence Nd^'^'^, Pm^''+ and Sm^^+ (cm^^). 



Ion 



Nd" 



Pmi''+ 



Sm^ 



60 



61 



62 



Level 



Energy 



5Sl/2 








4/5/2 


59279 


106000 


4/7/2 


64654 


111000 


5Sl/2 








4/5/2 


3973 


120000 


4/7/2 


10351 


126000 


4/5/2 








4/7/2 


7480 


7000 


5si/2 


56272 


-136000 



50 



55 



60 



65 



FIG. 2. Dirac-Fock ionisation energies of 5s (solid) and 4/7/2 
(dashed) levels for the Ag isoelectronic sequence. 



ionization energies of ~ 270 eV, the difference between 
them is very much smaller. In Table |T] we see that the 
excitation energy for ions near the crossing point is well 
within the range of optical lasers. 

Higher- Z ions with many valence electrons can have 
more complex behaviour as ion charge is increased. For 
example, in the neutral Th the 5s electron is a core elec- 
tron while the 5/ electron is a valence electron: therefore 
the 5/ orbital is above the 5s orbital on the energy scale. 
Moreover, it even lies above the 7s orbital. As seen from 
the experimental spectrum of the energy levels of the neu- 
tral (four-valence) Th, the energy of the 6(P7s5f state is 
above the energy of the 6(P7s^ state. In the hydrogen-like 
Th the energy of the 5/ orbital is approximately equal to 
the energy of the 5s orbital, i.e. it is significantly lower 
than the energy of the 7s state. Therefore, there should 
be an ion charge at which the 5/ orbital "crosses" the 
7s orbital, while at some higher charge it crosses the 6s 
orbital, etc. 



III. METHOD 

To calculate energy levels we use a combination of 
the configuration interaction and many-body perturba- 
tion theory methods (CI-hMBPT) ^ which has been 
extensively described in previous papers [26J, |27[ ; here we 
provide only a brief outline. Once we have level energies 
we can obtain a-sensitivity, by repeating the calculation 
for X = —0.01 and x = 0.01 {x is defined in (O). q is then 
extracted from the gradient of the transition energies ([6]) . 

The CI-t-MBPT method starts with the Dirac-Fock 
calculation for the closed-shell core. In all cases pre- 
sented here the core includes all shells up to Aspd. In the 
Ag isoelectronic sequence, this corresponds to the V^^^ 
approximation; for Sm^^+ it corresponds to V^'^'^. A 
single-particle i?-spline basis [2^ is then constructed in 
the potential of the nucleus and core electrons, which in- 
cludes valence orbitals and a large number of excited and 



virtual orbitals. In this paper we use a basis ISspdfg for 
the CI and SOspdfgh for the MBPT. 

For the two- valence-electron case, Sm^^"*", we perform 
the full CI calculation in the frozen-core approximation. 
Here the many-electron wavefunction is expressed as a 
linear combination of Slater determinants |/): 

^p = Y^ci\i) 

where the coefficients C/ are obtained from the eigen- 
value problem 
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HjjCj — ECi 
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and H is the CI Hamiltonian. Core-valence effects are 
taken into account using an MBPT operator S which is 
added to the CI Hamiltonian (see Section HI of [27i|). 
We calculate E to second order in the residual Coulomb 
operator, leading to the modified eigenvalue problem 
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The states \AI) include all Slater determinants of the 
single particle basis that have core excitations. Goldstone 
diagrams and analytical expressions for the one- valence- 
electron and two- valence-electron contributions to S are 
given in |27| . For the one- valence-electron ions presented, 
the CI calculation is unnecessary and we simply add the 
one- valence-electron part of the MBPT operator S to the 
Dirac-Fock energy. 



IV. RESULTS AND DISCUSSION 

In Table U we present the results of our method for 
ions of the Ag isoelectronic sequence near the 5s - 4/ 
orbital crossing point (see Section |ll]) . It is seen that 
even though the ionisation energies of these levels are 
very large (~ 270 eV), which is reflected in the very large 
q values, the transitions themselves can be within the 
optical regime. High sensitivity to variation of a can be 
achieved in clocks by comparing the very sensitive 4/ -^ 



TABLE II. Energy levels, sensitivity coefficients (q) (cm ^) 
and lifetimes (r) for lower states of Sm ^ . 



Configuration 



J 



Energy 



5s4/ -"F" 


2 












3 


1814 


987 


920 s 




4 


7282 


6350 


4.2 s 


Ap^H 


4 


495 


129898 






5 


6288 


135389 






6 


12067 


140213 




Af'F 


2 


9771 


131140 


16 ms 




3 


13746 


135303 


12 ms 




4 


14377 


134910 


35 ms 


5s4/ ^F° 


3 


13047 


6819 


50 ms 


Af'G 


4 


21850 


141015 


11 ms 


bs^'S 





28248 


-124689 





5s transition in Sm^^+ to the 4/ fine-structure transition, 
which in this case would be the anchor. Alternatively, 
if it is convenient to compare frequencies between ions 
with high accuracy, comparison of the 4/ — )■ 5s negative- 
shifting transition in Sm^^+ with one of the large positive 



shifting 5s — > 4/ transitions in Pm 



14+ 



Ndi3+ would 



allow an even higher sensitivity to a- variation. 

One of the most interesting cases is Sm^'*+. It has 
two valence electrons above closed shells which are in the 
5s and 4/ states for the low-lying configurations. The 
states of all three configurations 4/5s, 4/^ and 5s^ are 
relatively close to each other on the energy scale and are 
probably readily accessible to modern lasers. The transi- 
tions between these states correspond to the 4/ — >■ 5s or 
5s -^ 4/ single-electron transitions, which ensures strong 
sensitivity to the variation of the fine-structure constant. 
The results of calculations for Sm^'*'*' are presented in 
Table [TTl We are unaware of any experimental data for 
this ion, therefore only theoretical values are presented. 
Lifetimes are calculated using electric dipole transitions 
only (Table Uni). 

Consider, for example, the levels marked as GS 
(ground state), A (495 cm"!) and B (28248 cm^i). Both 
of these states have long lifetimes (no allowed dipole tran- 
sitions) and the linewidths of transitions from the ground 
state to A or B are probably sufficiently narrow to en- 
sure accuracy similar to that which is achieved in atomic 
clocks. If a varies in time, states A and B will move in 
opposite directions. Both of these states have some of 
the largest g-values ever seen in an atomic system. 

The relative change of the ratio of the frequencies of 
two transitions can be written as 



A{u}i/U}2) 
iuJl/uJ2) 



a 



where 



K = 



2q 



(10) 
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TABLE III. Reduced electric-dipole matrix elements for tran- 
sitions between lower states of Sm^^^ (a.u.) 



Level 1 
Configuration 



4f '-"F 

4/2 3p 
Af ^G 



5s4/ -^F" 
5s4/ ^F" 



5s4/ ^F° 



Level 2 
Configuration 
5s4/ -"F" 
5s4/ -^F" 
5s4/ '^F° 
5s4/ '■^F" 
5s4/ '■^F" 
5s4/ ^F° 
5s4/ '■^F" 
5s4/ ^F° 
5s4/ ^F" 
5s4/ '■^F" 
5s4/ '■^F" 
5s4/ ^F" 
4/2 '^H 

4/2 3j^ 

4/2 Sir 

4/2 3jc^ 
4/2 3^ 



Amplitude 

(2|MI|1) 



0.011800 
0.007480 
0.002030 
0.012300 
0.006650 
0.005530 
0.004450 
0.015500 
0.017900 
0.001850 
0.006820 
0.018100 
0.001280 
0.005890 
0.001970 
0.001840 
0.000078 



A (495 cm-i) and B (28248 cm^i): 



Ka ^ 525, 



Ku = -i 



A(u;a/^b) ^ 534^ 
{to A I lob) a 



(12) 



Substituting numbers from Table HIl one obtains for levels 



This is more than two orders-of-magnitude higher rela- 
tive sensitivity than that of the Hg+ frequency standard 
[K — —3.19), the system in which the current strongest 
constraint on the present-day time variation of a was 
obtained [9|. The very high sensitivity Ka comes from 
large q and small uj, and we should note that our theo- 
retical transition frequencies result from a cancelation of 
energy levels at the 0.01 - 0.1% level, and may only be 
accurate to perhaps ^ 2000 cm~^: enough to reduce the 
relative sensitivity Ka considerably. On the other hand 
the g-values in Table |lT] are very stable in our calcula- 
tion. Comparison of any of the lines in Table |lT] with B 
will give values of gi — (72 ~ 260 000 cm~^. This repre- 
sents an absolute enhancement in a-sensitivity five times 
that of the Hg+ZAl^ clock comparison. 

The lifetimes depend more strongly on the transition 
frequencies than on the dipole amplitudes. On the other 
hand, the accuracy of the transition frequencies may be 
low since they are obtained as the difference of the two- 
electron removal energies of two levels. Therefore, along 
with the lifetimes of the levels, we also present the calcu- 
lated values of the transition amplitudes (Table IIIip . If 
frequencies are measured then the lifetimes can be recal- 
culated using these amplitudes. Note that the amplitudes 
are small. This is because they correspond to s-f single- 
electron transitions, which cannot be an electric dipole 
transition. Therefore, the El transition amplitudes in 
the two-electron states are due to configuration mixing 
with suitable states, and this mixing is small. 



CONCLUSION 



We have shown that highly- charged ions present oppor- 
tunities for optical transitions with very high sensitivity 
to a-variation. The enhancement is proportional to the 
ionization energy of the states, which increases with ion 
charge as ^ {Zi + 1)^. We have presented a method to 
identify good experimental candidate ions that have suit- 
able transitions within the range of optical lasers. Apply- 
ing our method to the Ag isoelectronic sequence allowed 
us to identify Sm^''+ as having suitable transitions with 
the highest g-values seen to date in an atomic system. 
Furthermore it has a very high relative sensitivity, and is 



therefore an excellent candidate for studies of temporal 
a- variation. 

The ideas presented in this letter may be extended to 
other isoelectronic sequences, Au, Hg, etc., which would 
benefit from a larger Z^ enhancement, along with the 
{Zi + 1)^ enhancement we have discussed. 
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